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EFFECT OF STEADY VELOCITIES ON THE DYNAMIC RESPONSE OF 


LIQUID JET ATOMIZATION* 
by Marcus F. Heidmann 
Lewis Research Center 

SUMMARY 

The effect of steady axial, tangential, and radial velocities on the dynamics of liquid 
jet atomization in a traveling transverse acoustic mode is analyzed to examine the influ- 
ence of steady velocities on the stability of rocket engine combustors. The analysis as- 
sumes an atomization model prescribing jet breakup when jet distortion caused by aero- 
dynamic forces attains a critical value. 

An increase in the axial velocity difference between the gas and liquid jet degrades 
the response of atomization rate to pressure oscillations. The effect of lateral veloci- 
ties on response is more pronounced. Small tangential velocities in the direction of the 
acoustic wave travel and small radial velocities either toward the wall or toward the 
center significantly increase the magnitude of the response. Tangential velocities 
counter to the wave travel cause a substantial reduction in response. 

The results imply that variations in injection velocities or axial gas velocities can 
alter stability properties of rocket combustors. Stability properties should also be sen- 
sitive to small lateral velocities near the injector caused by nonuniform propellant in- 
jection or acoustic damping devices and by the acoustic streaming velocities associated 
with acoustic modes. 


INTRODUCTION 

In previous analyses (refs. 1 and 2), the dynamic response of liquid jet atomization 
to acoustic oscillations associated with traveling, standing and radial transverse acous- 
tic modes was presented. An atomization model (ref. 3) prescribing jet breakup when 
jet distortion caused by aerodynamic forces attains a critical value was used. The per- 

*Presented at Sixth ICRPG Combustion Conference, Sept. 9-11, 1969, Illinois 
Institute of Technology, Chicago, 111. 



turbating or acoustic particle velocities were found to affect significantly the in-phase 
response factor (mass release in-phase with the pressure perturbation) of the process. 
Particle velocities in-phase with the pressure perturbation increased the response factor 
above the pressure sensitive value, whereas out-of-phase velocities were suppressive. 
The analytical solutions obtained for these nonlinear (pressure amplitude dependent) ve- 
locity effects have provided some insight into the parameters which can contribute to the 
stability of rocket combustors. 

The analytical procedures used to extract particle velocity effects are also adaptable 
to examine steady velocities superimposed on the acoustic terms. Some steady velocity 
effects were previously reported. In the present study, a more general analysis of 
steady velocity effects on liquid jet response is presented. Axial, tangential, and radial 
velocities are treated. Principle extension of the previous analysis is with regard to 
steady radial velocities. 

Steady velocities are often suspected and sometimes known to affect the stability 
characteristics of rocket engine combustors. Instability has been initiated in research 
combustors by inducing a steady tangential velocity as in references 4 and 5. Tangential 
and radial gas injection from the wall has also been used to trigger instability in full- 
scale combustors (ref. 6). Such gas injection is usually most effective in the region 
near the injector, where propellants are atomized. In addition to these induced veloci- 
ties, steady flow effects on stability are often indirectly implied. Changes in stability 
behavior with injection velocity and contraction ratio may be examined from the view- 
point of changes in steady axial velocities. In some instances, recirculation near the 
injector is suspected to affect stability (ref. 7). Propellant distribution patterns are be- 
lieved to be responsible for recirculation and a variety of steady transverse velocities. 
Baffles, acoustic liners, and circumferential slots may also alter steady flow patterns 
and thereby influence stability characteristics. 

Because of the variety of steady flow conditions which can exist in a rocket combus- 
tor, this analysis attempts to provide some insight into how such steady flows can affect 
the dynamic response of the atomization process. The analysis has been directed toward 
steady velocities in traveling transverse acoustic modes. General analytical solutions 
are derived and specific evaluations are presented for some representative operating 
conditions. 


ANALYSIS 

The atomization model used is identical to that for the previous analyses of jet 
atomization (refs. 1 and 2). Aerodynamic pressure acting on an element of jet length 
causes an internal acceleration of liquid which eventually leads to breakup (see fig. 1) 
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If breakup is assumed to occur when displacement or distortion attains a critical value, 
the criterion for breakup of a continuous flow of jet elements is 
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f(t)dt** 


dt* = Constant 
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and the perturbation in atomization rate w, caused by fluctuations in aerodynamic load- 
ing as derived in reference 1, is 



dt Ti(t - r) 
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(All symbols are defined in the appendix. ) The parameter V is the magnitude of the 
relative velocity vector between the gas and the liquid. 

For this analysis, steady radial and tangential velocities are assumed to be super- 
imposed on the first order acoustic expressions for the traveling transverse mode. The 
modified acoustic properties are 
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In the notation used, u r is a positive steady radial velocity toward the wall, and 

u n is a positive tangential velocity in the direction of the acoustic wave rotation. 
u 2 

The expression for the p^V force defined by the acoustic properties of equation (4) 
can be reduced to a harmonic series. Disregarding second and higher harmonic terms, 
the expression is 
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Following the procedures of reference 2, the integration indicated in equation (2) is 
performed by assuming the difference between r and r to be negligible to give the 
atomization rate w’. The in-phase response factor N R (real part of w' with respect to 
P’) which is defined by 
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is evaluated to first order in P'. The solution for the in-phase response factor, includ- 
ing the effect of a distributed breakup process described in reference 2, is then 
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The first term within the brackets of equation (10) is identical to that obtained in the 
previous analysis of reference 2 where the in-phase response was given by 
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This converging periodic function with cot (eq. (12)) is shown in figure 2. It character- 
izes the response for conditions having perturbating particle velocities, steady axial ve- 
locity differences, and steady tangential velocities in the absence of steady radial veloc- 

A 

ities. Changes in steady and acoustic properties affect the value of (eq. (7)) and 
thereby affect the magnitude of the response. The form of the response function, how- 
ever, is not affected. 

A steady radial velocity introduces the second term within the brackets of equa- 
tion (10). The parameter (eq. (8)) is directly proportional to the steady radial ve- 
locity. The addition of this second term can appreciably alter the characteristic proper- 
ties of the response function. The extent of such property changes depends on the mag- 

A A 

nitude of Vg relative to Vj. Figure 3 shows the response properties for both positive 
and negative values of Vg for the ratio Vg/Vj. hi both cases, the response is en- 
hanced; that is, maximum and minimum are increased. The increase in the initial peak 
value is much larger, however, for positive values (radial flow toward the wall) than for 
negative values (radial flow toward the center). The increase in peak values is also ac- 
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companied by a shift in characteristic time c of at which these peaks occur. 

In the subsequent discussion, radial velocities which significantly effect the re- 
sponse function for typical operating conditions are presented. The quantitative effects 
of steady axial and tangential velocities on the in-phase response factor are also dis- 
cussed. 


DISCUSSION 

The analytical solutions for jet dynamics involve a multiplicity of interrelated pa- 
rameters, and a variety of interpretations relevant to rocket combustor stability could 
be extracted. The objective of the analysis, however, is to provide some insight into 
the dominant properties of jet dynamics rather than the intricate effects of an interplay 
of parameters. In pursuing this objective, precise modeling of jet atomization has been 
sacrificed to a degree that analytical solutions which retain dominant properties are ob- 
tained. 

The reliability of the analytical solutions decreases rapidly for characteristic times 
c of greater than 3tt or 4 tt. With large times the effect of reversing acoustic forces 
acting on the jet, restraining forces for liquid distortion, and the distributed properties 
of the breakup process must be considered in greater detail to give quantitatively signif- 
icant results. Such considerations will decrease the coherency of the jet mass and 
thereby reduce the response to insignificant levels as wr becomes large. 

The response in the region of off near n radians is relatively insensitive to the 
model limitations and retains a degree of quantitative significance. The discussion of 
steady velocity effects on the dynamics of jet atomization, therefore, will be confined to 
the initial peak value of the in-phase response factor which occurs near an off of n 
radians. This peak represents the maximum potential response for most jet conditions 
and, therefore, the magnitude of this peak implies the probability of unstable combustion 
being caused by atomization dynamics. 


Axial Velocity Difference 

The in-phase response factor for jet atomization (eq. (10)) is a function of the axial 
velocity difference between the liquid jet and surrounding gas. The peak value of the re- 
sponse factor in the absence of any steady radial and tangential velocities is very nearly 
given by the value of the function at off equal to 7 r radians: 
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Figure 4 shows the effect of axial velocity difference on peak response at a wall 
position (u r = 0). Increasing the velocity difference decreases the peak response. The 
effect is pressure amplitude dependent. Large velocity differences are required to 
cause suppression at high amplitudes. The effect changes with radial position is shown 
in figure 5 for the radii which divide the cross-sectional area of the chamber into four 
equal parts. A peak-to-peak pressure amplitude of 10 percent of mean pressure, 

P f = 0. 05, is used in the comparison (somewhat above the noise level in combustors and 
a lower limit for realistic mode amplitudes). As shown in figure 5, peak response is 
relatively unaffected by axial velocity difference near the center. An increase in axial 
velocity difference, however, is generally suppressive with regard to peak response and 
potentially more conducive to stable combustion. In effect, high axial velocity differ- 
ences tend to shield the jet from acoustic oscillations. 

Changes in design or operating conditions which would increase the axial velocity 
difference depend on the specific conditions being altered. For gas-liquid concentric- 
tube injection, the gas injection velocity can usually be made large relative to the liquid 
jet velocity and should provide potentially more stable combustion. This has been dem- 
onstrated in an extensive study with gaseous hydrogen and liquid oxygen combustors 
(ref. 8). High hydrogen velocities had a stabilizing affect which could override all other 
stability variables. For liquid- liquid injection, the liquid injection velocities often ex- 
ceed the prevailing gas velocity in the atomization zone. With such conditions an in- 
crease in the injection velocity or a decrease in the gas velocity by increasing contrac- 
tion ratio would tend toward higher velocity differences and less potential response. 

The stabilizing effect of decreasing contraction ratio was also observed in the hydrogen 
oxygen study (ref. 8). Frequently, however, flows near the injector are either acceler- 
ating or decellerating and it is difficult to formulate explicit statements about changes in 
axial velocity differences. Axial staging, where one propellant is atomized and vapor- 
ized more rapidly than the other, adds another variable to axial velocity control which 
may be useful in establishing large velocity differences for one propellant. 


Tangential Velocity 

A steady tangential velocity can either increase or decrease the peak response of jet 
atomization depending on the direction of the steady flow. In the absence of a steady 
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radial flow, the peak response is approximated by 
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Figure 6 shows the effect of both positive and negative tangential velocities on peak 
response at a wall position and for an axial velocity difference of 100 feet per second. 
Positive velocities (in the direction of wave rotation) increase the response, whereas 
negative velocities (counter to wave rotation) reduce the peak response. Steady veloc- 
ities as low as 5 feet per second can either double the response factor or reduce it to 
zero. At higher velocities, extremely large (compared with 1) positive and negative re- 
sponse factors for positive and negative velocities, respectively, are predicted. The ef- 
fect is similar at all radial positions as shown in figure 7 for the quadripartite area 
radii and a peak-to-peak pressure amplitude of 10 percent of mean pressure. A steady 
tangential velocity in the direction of the wave rotation, therefore, is potentially desta- 
bilizing and a velocity opposed to the wave rotation is stabilizing with regard to com- 
bustor dynamics. It should be noted that the large negative responses at an oof of n 
are accompanied by positive responses at an oof near 2ir. Small negative tangential 
velocities, however, will give zero response at all values of oof. 

Self-induced steady tangential velocities in rocket combustors have received little 
study. Some evidence of vorticity in solid propellant combustors has been observed 
(ref. 9), as predicted by "acoustic streaming" considerations in acoustic theory. Ref- 
erence 10, for example, predicts that low order traveling modes should induce a steady 
wheel flow in a direction opposite to the wave rotation. This would have a stabilizing ef- 
fect on the atomization process. Such wheel flow will not develop in regions of high 
Mach number axial flow and little emphasis has been placed on such flow in liquid pro- 
pellant combustors. These "acoustic streaming" or steady tangential velocities are re- 
lated to the boundary layer conditions at the wall. Unusual boundary layer conditions 
created by circumferential slots or acoustic liners in the wall may be expected to affect 
"acoustic streaming" - especially near the injector where axial Mach number is low or 
negligible. The development of small steady tangential flows and the effect of slots and 
liners on such flows should not be neglected considering the extreme sensitivity of jet 
dynamics to such flows. 

The effect of steady tangential velocities (induced by tangential injection of a sec- 
ondary gas) on jet dynamics and combustion stability has been demonstrated in research 
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combustors (refs. 4 and 5) and analyzed in reference 1. Introducing a tangential velocity 
by secondary injection, canting of the injection pattern or other means, into an atomiza- 
tion dependent combustor will usually have a destabilizing effect, independent of direc- 
tion. In a developing acoustic mode, the preferred direction of wave travel is in the di- 
rection of giving maximum response. Prolonged steady flow, therefore, will cause a 
wave motion in the direction of the steady flow. Momentary flow opposed to an estab- 
lished wave, however, may be useful in suppressing instability as discussed in refer- 
ence 11. 


Radial Velocity 

As discussed previously and shown in figure 3, a steady radial velocity either to- 
ward the wall or toward the center can increase the peak response. An analytical solu- 
tion at a fixed value of cot cannot be used to express peak values as was the case for 
axial and tangential velocities. Radial velocities significantly affect the value of c or at 
which the peaks occur and direct analytical solutions for peak values must be used to ex- 
amine the response for typical flow conditions. 

Figure 8 shows the effect of both positive and negative radial velocities on peak 
response at a radius of 0.71 (area median radius) and an axial velocity difference of 
100 feet per second. Figure 9 shows the effect with radial position. The amplifying ef- 
fect of radial velocities is most severe at low amplitudes and near the center of the 
chamber. Flow toward the wall also has a much larger effect than flow toward the cen- 
ter. Flow toward the wall has an effect comparable to that for positive tangential veloc- 
ities, whereas the increase in response is reduced to about one-third these values for 
flow toward the center. For both flow directions, however, the increase in peak re- 
sponse in substantial and small radial velocities could significantly alter the stability 
characteristics of rocket combustors. 

Steady radial velocities are not predicted from acoustic theory for "hard wall" 
cylinders as was the case for tangential velocities. Slots and linear, however, provide 
"soft wall" conditions which may allow radial velocities to develop - especially for the 
zero or negligible Mach number flow conditions at the injector face plate. The nonlinear 
acoustic solutions for such flows have received little attention. 

Steady radial velocities can develop in actual rocket combustors from other causes. 
Recirculation of combustion gases can introduce steady radial flows sweeping across the 
injector face. An axial-radial recirculation path may be expected whenever a void exists 
between the injection pattern and the chamber wall. The direction of this flow may be 
uncertain. Radial flow toward the wall may be expected if combustion occurs close to 
the injector face (see fig. 10). The opposite flow direction is possible when combustion 
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occurs beyond some critical distance from the injector. Such flow patterns could signif- 
icantly affect jet response and combustor stability but have received only qualitative at- 
tentions as in reference 7. 

The influx and discharge of gases from slots and acoustic liners may also be in- 
strumental in creating or affecting radial flow patterns. For example, the flow process 
for acoustic liners is usually considered to be nonlinear for large pressure amplitudes, 
as shown in figure 11. Potential flow occurs during influx, but the flow is dissipative 
with respect to kinetic energy during discharge. The spacial differences in the two flow 
patterns can cause recirculation of combustion gases during periodic flow and directly 
affect jet dynamics. This pumping action of acoustic cavities can also modify existing 
radial flows related to nonuniform mass distribution and other causes and thereby indi- 
rectly affect jet behavior. The sensitivity of the atomization process to small radial 
flows places increased emphasis on design features which affect these previously ne- 
glected flow properties. 


Nonlinear Instability 

The dynamic response of the jet atomization process varies with pressure ampli- 
tude and, therefore, provides some insight into the triggering of nonlinear instability in 
rocket combustors. This was discussed in the previous analysis (ref. 1), but the effect 
of steady velocities on dynamics provides some additional concepts on nonlinear behavior. 

Previously, it was argued that a pressure disturbance reduced the amplitude sensi- 
tive characteristic time air, as expressed by equation (11). A disturbance which is suf- 
ficiently large to reduce wr to a value less than n ( see fig. 2) could trigger instability 
during the decay of the input disturbance. During decay, c of increases. In-phase re- 
sponse also increases and instability would develop if and when the in-phase response 
exceeds the acoustic losses of the system. 

An input disturbance (bomb or shock with mass addition) which also introduces a 
steady radial or tangential velocity modifies the concept of triggerable instability. De- 
pending on the direction and magnitude of these steady flows, they may either increase 
or decrease the probability of sustained instability. A transient flow which enhances the 
peak in-phase response (radial or in the direction of wave travel) may persist sufficiently 
long to allow the disturbance to grow rapidly into self-sustained instability by providing a 
response magnitude which drastically exceeds acoustic losses. Similarly, the probabil- 
ity of instability would be decreased by transient flows which suppress the response. 

Circumferential slots and acoustic liners may also generate momentary flow condi- 
tions and affect the probability of triggering nonlinear instability. Of particular concern 
are wall cavities located near the injector. These cavities can be charged with inert or 
combustible gases by an input disturbance and discharge during the decay period. Both 
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the influx and discharge could create momentary lateral velocities in the low Mach num- 
ber flow regime near the injector and affect the stability rating. This behavior may 
contribute to the stabilizing affects obtained from slots and liners during stability testing. 


CONCLUDING REMARKS 

This analysis of jet atomization shows that steady low-level radial and tangential 
velocities could significantly affect the stability of rocket engine combustors. The sig- 
nificance of steady velocity on stability, however, is not limited to jet atomization. 

Most atomization processes (impinging jets, spray nozzles, etc.) are dependent on aero- 
dynamic factors in a manner similar to that for jet atomization. Although the analyses 
for other atomization methods is more complex than for jets, the dependence on aerody- 
namic forces would give a dynamic response similar to that for jets. 

Accepting the validity of the jet analysis and its qualitative adaptation to most injec- 
tion methods places increased importance on low-level lateral gas flow in the region of 
the injector. Some sources of lateral flow were postulated and discussed in this report. 
These and others must be more fully explored with regard to the magnitude, and the pa- 
rameters affecting lateral flow before the relevance of steady flow effects to actual com- 
bustors can be established. 

Some studies of lateral flow behavior which could influence the design concepts for 
stable rocket combustors are: 

(1) "Acoustic streaming" velocities for "soft, " dissipative and irregular walls with 

emphasis on boundary layer variables that enhance tangential and radial veloci- 
ties which suppress atomization response 

(2) Lateral and axial mass distribution effects on lateral and recirculation velocities 

with emphasis on near-uniform mass distribution rather than isolated sources 

(3) Influence of acoustic liners, slots, and cavities on lateral and recirculation ve- 

locities caused by mass distribution effects under dynamic conditions 

(4) Transient lateral velocities caused by slow discharge of abruptly charged wall 

cavities with attention on the probability and consequences of charging with com- 
bustible mixtures 

(5) Lateral velocities for combustion chambers with minor deviations from cylindri- 

cal geometry such as concave and convex injector faces, tapered, and contoured 
chamber walls 

(6) Steady flow patterns in combustion chambers with baffles and surface protrusions 
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SUMMARY OF RESULTS 


The analysis of steady velocities effects on the the dynamic response of jet atomiza- 
tion in traveling transverse acoustic modes has shown the following with regard to the 
peak in-phase response to pressure oscillations: 

1. High axial velocity differences between the liquid and surrounding gas suppress 
the response. 

2. A steady tangential velocity in the direction of the acoustic wave rotation in- 
creases the response, and a counterwave velocity suppresses the response. * 

3. A steady radial velocity increases the response with flow toward the wall giving 

about three times the increase obtained from flow toward the center. t 

4. The sensitivity of response to tangential and radial velocities is sufficiently large 
to be of concern with regard to velocities caused by nonuniform propellant distribution, 

’’acoustic streaming, ” and acoustic damping devices in actual combustors. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 25, 1969, 

128-31. 
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APPENDIX - SYMBOLS 


speed of sound, ft/sec 

jet diameter, ft 

Bessel function of first kind of 
order n 

response factor, real or in- 
phase component of w'/P’, 
dimensionless 
2 

pressure, lb/ft 

radial position, dimensionless 

time, sec 

velocity, ft/sec 

magnitude of relative velocity 
vector, ft/sec 

parameters defined by eqs. (7) 
and (8) 

mass flow rate, lbm/sec 

radius parameter, 1. 841 r 

ratio of specific heats 

critical value of relative 
distortion, 5/D 


6 jet distortion, ft 

(jo frequency, rad/sec 

0 phase angle, rad 

p density, lb/ft^ 

a instantaneous deviation of 

from mean 

r jet breakup time, sec 

Subscripts: 
g gas 

1 liquid 

r radial component 

rms root mean square 

z axial component 

6 angular component 

w wall 

Superscripts: 

( — ) mean or average value 

(") maximum value 

( )' perturbation about mean 
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Figure 1. - Jet atomization model. 
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Maximum in-phase response, tNR> max Maximum in-phase response, %) max 



Figure 5. - Steady axial velocity-difference effect on response with 
radial position. Pressure amplitude, 0.05; speed of sound, 5000 
feet per second; ratio of specific heats, 1.2. 



Figure 6. - Steady tangential velocity effect on response with 
pressure amplitude. Radius, 1.0; |uj - u z |, 100 feet per 
second; speed of sound, 5000 feet per second; ratio of 
specific heats, 1.2. 
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Figure 7. - Steady tangential velocity effect on response 
with radial position. Pressure amplitude, 0.05; 

|u t - u z |, 100 feet per second; speed of sound, 5000 
feet per second; ratio of specific heats, 1. 2. 


Pressure 

amplitude, 

P' 



Figure 8. - Effect of radial velocity on response with ampli- 
tude. Radius, 0.71 (area median); |u^ - u z |, 100 feet per 
second; speed of sound, 5000 feet per second; ratio of 
specific heats, 1.2. 
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Figure 9. - Effect of radial velocity on response with radial 
position. APp_p , 10 percent; Iuj - u z |, 100 feet per 
second; speed of sound, 5000 feet per second; ratio of 
specific heats, 1.2. 
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Figure 10. - Radial winds 
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mass distribution effects. 
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Figure 11. - Recirculation from nonlinear cavity flow processes. 
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